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A localized spin coupled to the spins of a conduction 
electron system via a Heisenberg exchange interaction 
is known to give rise to the Kondo effect, provided the 
coupling is antiferromagnetio^. At sufficiently high tem- 
peratures or energies, well above the characteristic en- 
ergy scale referred to as the Kondo temperature, Tr-, the 
signature of the Kondo effect is a logarithmic variation 
of various observables with temperature (or other ener- 
gies) . As first demonstrated by Kondo^ , such logarithmic 
behavior appears already in perturbation theory to low 
order in the exchange coupling. At low temperatures, 
however, for T <^ Tk, the local spin is screened by the 
conduction electron spins, and the system enters a local 
Fermi liquid state, characterized by integer temperature 
power laws. 

The many-body resonance state forming near Tk is 
comprised by infinitely many virtual particle hole excita- 
tions leading to a peak in the conduction electron scatter- 
ing amplitude at the Fermi energy: the Kondo resonance. 
In equihbrium systems the Bethe ansatz method allows 
to analytically calculate thermodynamic properties'^, and 
dynamical properties can be determined with the help of 
Wilson's numerical renormalization group^. These meth- 
ods reveal the universal nature of the Kondo effect: a sin- 
gle characteristic scale Tk determines the physics and all 
physical quantities are universal functions of, e.g., T /Tk- 
Already in 1970 Anderson^ suggested a simple and effi- 
cient method, known as "poor man's scaling" or pertur- 
bative renormalization group, to resum the leading loga- 
rithmic terms in perturbation theory and to establish the 
scaling behavior. It requires nothing more than low order 
perturbation theory, and provides a controlled approxi- 
mation for T > Tk, i.e. as long as the running exchange 
couplings remain small. 

While the Kondo effect was discovered in metals con- 
taining magnetic impurities in the 1960's (with experi- 
mental observations dating back to the 1930's), and most 
of the theoretical development took place before the mid 
1980's, it had a revival in the 1990's in the context of elec- 
tron transport through quantum dots weakly coupled to 
leads. Provided that the dot carries a net spin, a Kondo 



resonance develops, thus admitting resonant tunneling 
of electrons^. This leads to a removal of the Coulomb 
blockade, i.e. an increase of the conductance from small 
values up to the quantum limit, as the Kondo resonance 
develops. We shall refer to this type of quantum dot as 
a Kondo dot. It has been observed that the Kondo effect 
is quenched by raising the transport bias voltage V well 
above Tk, i.e. eV ^ Tk, and that the presence of a 
magnetic field splits the zero-bias conductance peak into 
two distinct peaks, located at bias voltages roughly equal 
to plus and minus the Zeeman splitting of the spin on the 
jQ^j8j9ji0jiiji2 Purthermore, it has proven possible to col- 
lapse the data from different Coulomb blockade valleys 
onto one universal conductance curveii*i2ii^, thus estab- 
lishing the aforementioned one-parameter scaling out of 
equilibrium. 

The observation of Kondoesque conductance anoma- 
lies, exhibiting a logarithmic temperature dependence of 
the zero-bias peak height and a Zeeman splitting of the 
peak in finite magnetic field, actually has a much longer 
history; it has often been observed in more traditional 
tunnel junctions involving tunneling via magnetic im- 
purities. In metal-insulator-metal junctions the Kondo- 
anomaly can result from magnetic impurities present in 
the insulating metal-oxide barrier, from surface states at 
the metal/metal-oxide interfac9i^ii^*iS4iLi& or even from 
unpaired electrons residing on organic radicals in a poly- 
merized benzene barrier—. The same type of conduc- 
tance anomaly has been detected in semiconductor-metal 
(Schottky) junctions^2iSi, in which the neutral shallow 
donors in the semiconductor depletion layer provide the 
spin-1/2 moments which incite the Kondo correlations. 

In 1966, AppelbaumSSiS^ and Anderson^i demon- 
strated that this type of conductance anomaly can result 
from so called exchange tunneling processes, in which an 
electron tunnels from one electrode to the other via an 
intermediate magnetic impurity orbital and at the same 
time flips its spin. This mechanism was shown to lead 
to a Kondo effect enhancing the charge transfer across 
the dot, and it explains in a natural way why the con- 
ductance is peaked at a bias voltage corresponding to 
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the Zeeman splitting of the impurity moment: the fi- 
nite bias has to supply the energy to flip the spin in the 
presence of a magnetic field. Appelbaum calculated the 
tunneling conductance from a simple golden rule expres- 
sion, in which the tunneling amplitudes were determined 
from third order perturbation theory, including the lead- 
ing logarithmic corrections reflecting the Kondo effect. 
His result did capture the qualitative features of the con- 
ductance anomaly, but any quantitative agreement with 
experiment has been restricted to the case of zero mag- 
netic field. Wolf and Losee^^ later suggested an exten- 
sion of Appelbaum's theory in which the logarithmic en- 
hancement was cut off by the lifetime broadening of the 
Zeeman split spin-levels. Indeed this extra feature im- 
proved the agreement with experiments somewhat, but 
even with several parameters, this did not allow for a fit 
to experiments at finite magnetic fielciiSiiLi2iS242&. 

The mechanism of exchange tunneling, suggested by 
Appelbaum and Anderson, was originally built into a 
Hamiltonian with constant, a priori unknown, exchange 
couplings to the spin. While this approach relied on an 
exchange tunneling term connecting the two electrodes 
(see the model (QJ below), a parallel development by 
Solyom and Zawadowski^^ suggested that the conduc- 
tance anomaly should arise from the energy dependent 
renormalization of the local density of states, induced 
by Kondo spins near one of the electrodes not directly 
coupled to the other side. Initially the two proposals dis- 
agreed even in the sign of the change in conductance, 
but these differences were resolved in a later work by 
Appelbaum and BrinkmanSS, using Zawadowskis alter- 
native approach to tunnelingS^iffi. The final reconcilia- 
tion of these ideas came with the work by Ivezic"^^, tak- 
ing the nonequilibrium Keldysh approach to tunneling, 
developed earlier by Caroli et alMi. The fact that the 
electrodes were out of mutual equilibrium was taken into 
account, and it was demonstrated that the earlier (essen- 
tially equilibrium) treatments^i^ were only correct in 
cases when the impurity was located nearby one of the 
electrodes, whereas an impurity somewhere in the middle 
of the barrier constitutes a true nonequilibrium problem. 
A good review of these earlier works may be found in 
Ref.H 

With the discovery of the Kondoesque tunneling 
anomaly in various types of quantum dots, this vener- 
able problem has recently been revived. Meanwhile, ex- 
perimenting with quantum dots offers much better con- 
trol over the parameters defining the problem, and this 
in turn allows for a more systematic study of the physics 
underlying the conductance anomaly. Since the zero-bias 
anomaly arises only in Coulomb blockade valleys corre- 
sponding to an odd number of electrons occupying the 
dot, the effective local moment can be ascribed to a single 
electron in the uppermost energy level. In this way the 
'magnetic impurity' spreads over the entire dot and there 
is no confusion as to whether the impurity is located close 
to, or even residing in, one of the electrodes or whether 
one should average over many impurities, issues which 



were all very important in the conventional tunnel junc- 
tions mentioned above. In particular, there is no reason 
to believe that the dot-spin should be equilibrated with 
one particular lead, which was pointed out by Ivezic— 
to constitute a true nonequilibrium problem. Within 
the Anderson model, this problem has been studied us- 
ing an equations-of-motion technique combined with the 
non-crossing approximation'^^ '^'"'"^'' ''^"^* ''^ and by means 
of nonequilibrium perturbation theory, expanding in the 
hybridization strength^. Other works have taken the 
Kondo model as their starting point, but have mostly fo- 
cused on the effects of an applied ac-bias in the case of 
zero magnetic field^ii^. 

A complete theory of the Kondo dot in a nonequilib- 
rium stationary state, i.e. in the presence of a finite cur- 
rent flowing through the dot, does not exist yet. Most of 
the methods which have proven so successful in dealing 
with the equilibrium Kondo problem, appear to have no 
trivial extension to a nonequilibrium situation. Schiller 
et alm^ have successfully applied bosonization techniques, 
originally devised by Emery and Kivelson-^ in their so- 
lution of the two-channel Kondo problem, to calculate a 
number of observables near a certain Toulouse point in 
the presence of a finite bias voltage. It is unclear, how- 
ever, to what extent these results apply to the generic 
Kondo model. Konik, Saleur and Ludwig^ started from 
the Bethe Ansatz solution in equilibrium to construct 
approximate scattering states, in the presence of a fi- 
nite voltage. To what extent their approximations and 
boundary conditions are valid is not obvious to us. To 
our knowledge, even the perturbation theory in a true 
nonequilibrium situation has not yet been worked out to 
leading logarithmic order in the presence of a magnetic 
field. This latter task will be taken on in the present 
paper, which presents a detailed analysis of the nonequi- 
librium perturbation theory to leading logarithmic order, 
generalizing Appelbaum's result to the case where the 
spin is not equilibrated with one of the leads. 

Our aim here is to understand the physical processes 
governing the nonequilibrium situation and to formu- 
late a starting point for the resummation of the lead- 
ing logarithmic terms in perturbation theory. Even in 
the perturbative regime, when magnetic fields, voltages 
or temperatures are large compared to Tk, such a re- 
summation is necessary not only to recover the correct 
universal scaling behavior but also to be able to fit exper- 
iments quantitatively^. In some of the early worksSiii, 
already the perturbation theory was resummed using 
approximations of Nagaoka"*^, Abrikosov"'^ or Suhl and 
Wong^", mainly as an attempt to capture the strong 
coupling regime. As for the quantum dots, Kaminski 
et alM. first suggested a poor mans scaling method to 
deal with this problem at finite voltage, but this ap- 
proach did not encompass the case of a finite magnetic 
field and therefore avoided the conceptual problem con- 
nected to the calculation of the magnetization. Recently, 
the present authors suggested a different poor man's scal- 
ing approach"*^, which effectively resums the logarithmic 
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corrections and recovers one-parameter scaling, even in 
the presence of a magnetic field. As an important check 
on this approach, it was crucial that an expansion in bare 
parameters would correctly reproduce the leading loga- 
rithmic corrections found from perturbation theory, and 
indeed this was found to be the case. Main results of 
the nonequilibrium perturbation theory presented here 
were briefly stated in Ref. [i^ The details of our scaling 
approach will be provided in a subsequent publication^. 

The remaining parts of the paper are organized as fol- 
lows. In section we define our model and briefly re- 
view the pseudo fermion formalism to be used through- 
out the paper. Section Hll reviews the diagrammatic rules 
of nonequilibrium perturbation theory using the Keldysh 
formalism, and introduces a few basic building blocks 
from which all subsequent diagrams will be constructed. 
In section IIIII we solve a quantum Boltzmann equation 
to find the nonequilibrium distribution function for the 
local spin, and demonstrate that this leads to rather dra- 
matic effects of the bias-voltage on the local magnetiza- 
tion. Finally, in section lTVI we derive an expression for the 
tunneling current which takes the nonequilibrium mag- 
netization as its input. Section Ivl contains a summary 
and discussion of our results. 



I. THE MODEL 

We consider the Hamiltonian 



(1) 



where hl,r = ieV/2, Sz is the spin-1/2 on the dot and r 
the Pauli-matrices. We shall use the dimensionless cou- 
pling constants gi = N{0)Ji, where N{0) is the density 
of states per spin for the conduction electrons and Jll, 
Jrr and Jlr = Jrl the real- valued exchange constants. 
For simple quantum dots in the Kondo regime, which can 
be described by an Anderson model, the exchange cou- 
pling constants are related by = gLLdRR- However, 
in more complex situations, like in double dot systems 
in the Kondo regime, no such relations exist and we will 
therefore treat as an independent parameter of our 
Hamiltonian Q. For notational convenience we employ 
the shorthand notation g^ = {gLL + gRR)/'^, and unless 
specifically stated otherwise, we will henceforth use units 
where h = ks = gUB = e = 1. 

In order to proceed with a perturbative calculation, 
which includes also the effects of a magnetic field, it is 
convenient to apply a fermionic representation of the lo- 
cal spin operators. We choose here Abrikosovs'*^ pseudo 
fermion representation, which, in terms of fermionic 
charge-neutral spin 1/2 operators /, reads 



j'^'-n'f-f' 



(2) 



Since only the singly occupied fermion states have any 
physical relevance, this representation must be supple- 
mented by a projection onto this physical part of the 
Hilbert space, effectively excluding doubly occupied and 
empty states. To this end, the pseudo fermion is en- 
dowed with a chemical potential A which is kept finite 
throughout the calculation, using a grand-canonical en- 
semble average. As demonstrated, for equilibrium and 
nonequilibrium situations, in Refs. and fssl respec- 
tively, the physically relevant, i.e. canonical ensemble 
averaged, expectation value of an observable O is ob- 
tained as the limiting value 



{0)q=i 



lim 



(Q) 



A 



(3) 



where Q = f^f-y is the pseudo fermion number opera- 
tor. Note that in the common case, where the observable 
O has zero expectation value in the Q = ensemble, one 
may leave out the Q-operator from the numerator in ||2J). 

Since A enters as a chemical potential, thermal av- 
erages taken with a finite A contain various powers of 
g-A/T_ particular {Q)x ^ e^-^^"^ and therefore the 
limit of A — !■ oo in Eq. J^J effectively picks out the terms 
in (OQ)x which are also proportional to e~^^^ . Likewise, 
in calculating any thermal average of interest at finite A, 
one is allowed to retain only the terms of lowest order in 



II. KELDYSH DIAGRAMMATICS 

In setting up the nonequilibrium perturbation theory, 
we shall comply with the conventions in Ref. All 
Keldysh-space matrix propagators arc represented in the 
usual upper-right triangular form 



G = 



G^ 



(4) 



and the individual entries will be denoted by Latin in- 
dices. From this basic Green function, one may obtain 
the usual spectral, lesser and greater functions as 



A 
G< 
G> 
Rc[G] 



(G^-G^ + G^)/2, 
(G^ + G^-G^)/2, 
{G" + G^)/2, 



(5) 
(6) 
(7) 
(8) 



77' 



implying the R superscript in the real part. A corre- 
sponding notation will be used for self-energies, except 
that the spectral function will be replaced by the broad- 
ening F = i(S]^ — T,^). We shall henceforth denote con- 
duction electron (ce), and pseudo fermion (p/) correla- 
tion functions by capital Latin, and calligraphic letters, 
respectively. 

Perturbation theory in terms of Keldysh matrix prop- 
agators involves the bare four-point vertex and the mea- 
surement vertices depicted in Fig. ^ These have the 
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FIG. 1: (a) The bare conduction electron-pseudo fermion 
vertex, (b) The absorption, and (c) the emission vertices. 
Propagators are dashed (p/) and full (ce) lines. Note that 
the absorption, and emission vertices could equally well have 
been drawn with p/-propagators. 
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FIG. 2: The conduction electron-pseudo fermion bubbles in 
Peterls (a) and Cooper (b) channels. Conduction electron 
propagators have lead-index a, and pseudo fermion propaga- 
tors carry the spin index 7. Latin indices refer to the Keldysh- 
matrix structure in ijl}. 
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Diagrams should be interpreted with the usual zero- 
temperature Feynman rules, including a prefactor 
(« J)"(— 1)^^^+^^" to any diagram of order J", having 
Fp/,ce closed pseudo fermion, or conduction electron 
loops. Both flavors propagate with spin-indices 7, cr = 
±1, to which we apply the convention that 7 = —7. Fur- 
thermore, the conduction electron Green functions carry 
a lead-index a € {L,K\. The Einstein summation con- 
vention is implied for all indices except where otherwise 
stated. 

Since we assume the total system to have reached a 
steady state, all Green functions depend only on one fre- 
quency. The bare pf propagator has the spectral function 



A^{uj) = 2tt5{uj + -iB/2). 



(12) 



and corresponding retarded and advanced Green func- 
tions. For the Keldysh component, we make the usual 
ansatz that 



(13) 



which defines the p/-distribution-fimction n~f\{Lij). At 
zero bias-voltage, where the spin is equilibrated with the 
conduction electrons, n^\{oj) reduces to the Fermi func- 
tion f{uj + \y= l/(e("+^'/^ -I- 1). We follow the conven- 
tion of Ref. |5J and place the chemical potentials in the 
distribution functions rather than in the spectral func- 
tions, which prevents the parameter A from pervading 
the formulas. The unprojected mean occupation num- 
bers are given by the integral 



' cLu 



(14) 



that is n^x = n^\{—^B/2) for bare A^. After taking the 
limit A — > 00 for a given observable and thus performing 



the projection onto the physical spin-states, one can use 
the constraint + = 1 to write the local magneti- 
zation as M — 2n| — 1. At times, we shall also use the 
shorthand M^^\ — 2n^^\ — 1 prior to projection. 

The bare conduction electron Green functions depend 
on momentum as well as frequency, but since the inter- 
action and the pseudo fermions are local in space, all 
internal lines in a Feynman diagram involve only prop- 
agation in time, and we are allowed to work with local, 
momentum integrated, ce-Green functions. Assuming a 
constant density of states in a band of width 21?, centered 
at zero, the momentum integrated spectral function reads 



A(cj) = 2-KN(Si)Q{D- |cj|). 



(15) 



with iV(0) = 1/(2_D). Assuming the electrons in separate 
leads to be in thermal equilibrium, the lead-dependent 
Keldysh Green function takes the form 



(w) = -iA(w) tanh 



2T 



(16) 



in terms of the chemical potentials fJ.L = ^IJ-r — V/2. 
The spectral function (|15|l . implies a small real part 

Re[G^(w)] = ^ln(^) w NiO)uj/D which we can 

safely neglect for most of our discussion. 



Second Order Vertex Functions 



The central objects in the diagrammatics of this problem 
are the ce — pf bubbles depicted in Fig. |21 insofar as 
they carry the logarithmic integrals which eventually lead 
to the Kondo effect. The bubbles in the Peierls- and 
Cooper-channels are evaluated as 

'^.Ktie) = J^Gt^{u; + e)g^;ic.) (17) 
"Etis) = I^Gi^{e-u;)gt;{u;), (18) 
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from which one can readily determine the various 
Keldysh-components by straightforward integrations: 



organize all entries in the Peierls channel as 





= ^K.A{e) = 0, 
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= iG^{e--fBl2), 
= -iG^{e-^B/2), 
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— ill 
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— ±.22 
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= -^Gf(e-7S/2) 




t" 

i22 


-7-22 
— ill 
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— il2 


- 721 

— i21 


= 0, 



-i7V(0)ln 



-Gf(£-7i?/2) 



AN{Q)\tl 



7^/2)2 +T2 



for the Peierls-channel, from which the bubbles in the 
Cooper-channel are obtained as 



7 



^ 



A,R 



(20) 



A,R,K 
K 



7^_R,a(£) 



In both channels the Kondo effect derives from the imag- 
inary part of bubbles with a if-component conduction 
electron together with a retarded, or an advanced pseudo 
fermion. In such combinations, Re[C/-y(a;)] oc l/w is con- 
voluted with tanh(e -I- ll)/2T) to produce the logarithmic 
enhancement. 

To work out higher order Feynman diagrams, it is con- 
venient to construct the second-order renormalized ver- 
tex functions, obtained by attaching a Keldysh vertex @ 
to each end of these bubbles: 



Ac' d a y^d'c' / ^\ \ cd' 
ab' l!Sdj'a'{£)J^a'b^ 

Acd' OtlPd C { ^\ AC d 
ab' l!Sdj'a'{£)J^a'b- 



cd' 



(21) 
(22) 



Working out the contraction of Keldysh indices one may 



(23) 



satisfying the symmetry 2^ = ~ J'^, and similarly 
in the Cooper-channel 



i22 


— il2 


— i21 


— ill 


. It-. 


T" 

il2 


— ill 


— i21 


— i22 


. It-. 


I" 
i22 


— ill 


— i21 


— il2 


4i 


I" 
ill 


— i22 


— il2 


— i21 


= 0, 



(24) 



-cd 



where X^jj = X^-;, — X^j. For clarity, we have temporarily 
suppressed the variables 7, a and e, and written only the 
Keldysh indices. Furthermore we have introduced new 

functions, X^'^'^ and X ' ', with 



jRA 
■R,A 



X 



rK 



i (jC 



R,A 
K 

R,A 



^R,A 



AM 
K 



I^A + ^r) 



T^ 



+ ^R + 



(25) 



satisfying the relation: ^l"^'^''^ (e) = -^I^'^'^ {-e). 

From their definitions, these new functions X^'^'^ are 
seen to have poles in the lower (R), the upper {A), or both 
(K) half-planes, and to satisfy that (X^(e))* = X^(£). To 
be specific, one finds from (|19|l : 



^X«(£) = -M^^xG^ie - 7S/2) - ^G^is - jB/2) + N{0) In , (26) 

^I^ie) ^ ^A(£-7S/2)|l-f Af,,Atanh(^^^i^^^^^)|. (27) 
The real, and the imaginary part of X^ are 

Re [^I{e)] - NiO) In ( _ _ ^^/g)^ + ) " ^'^^^^^'^ [^"^^ ~ ^^/^^^ ' ^^^^ 

Im [^1(8]] = iA(£ - 7B/2) (^Af^.A + tanh ~ , (29) 
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while the lesser and greater components take the following form: 



^X<(£) = 2iA{e - jB/2){l- n^^x)f{e - fic. - jB/2), (30) 

^7 V 



"I>{e) = 2iA{e~jB/2)n^^^{l-f{e-fi^-jB/2)). (31) 



Note that since Re[G(cj)] « N{0)uj/D for w < D, 
contributions from the last term in (|28() . of order 
max[B, V]/D <C 1, can safely be neglected. 

The real and imaginary parts (|28l I29() satisfy the 
Kramers-Kronig relation, and altogether their analyti- 
cal properties allow us to interpret X as genuine Keldysh 
Green functions describing the (time-) parallel or anti- 
parallel propagation of conduction electrons and pseudo 
fermions. The different components {R, A, <,>) of this 
mixed bubble could have been written down immediately 
using the Langreth rules--^ for analytical continuation, 
and in fact this method is also very convenient for de- 
termining the different components of the second order 
self-energy. For the third order self-energy, however, we 
find the Keldysh matrix structure to be more convenient 
when dealing with the large number of contractions ap- 
pearing to this order. The catalogues H23|) and H24|l have 
been worked out for this purpose and are used repeat- 
edly in working out the different Keldysh contractions 
appearing throughout the paper. 



III. NONEQUILIBRIUM MAGNETIZATION 
A. Pseudo Fermion Self-Energy 

With these basic diagrammatic objects at hand, we can 
now readily evaluate the Feynman-diagrams for the pf 
self-energies shown in Fig. O It is, however, essential to 
realize that the occupation functions n^x on the dot are 
completely undetermined to zeroth order of perturbation 
theory, i.e. in the absence of a coupling to the leads^. 
These are therefore kept as free parameters to be deter- 
mined later. The first order p/ self-energy vanishes, since 
we have not included a Zeeman term for the conduction 
electrons. Nevertheless, the corresponding Hartree term 
would be entirely real and contribute only by a constant 
shift of the pf energy levels (cf . appendix) . 



Second order self- energy 



Neglecting the Zeeman term for the conduction elec- 
trons, the a6-component of the second-order self-energy 
from Fig. O is given as 



where a summation over ce-spin has led to the tensor 



r;,,r;^,r;^,,r^,^=<5^y +2r^\. (33) 



Contracting the Keldysh indices c and d, using (|23|l , and 
using the analytical properties of X, the K component of 
this self-energy may be written as 

<?^(-) = ^-^JLj^mHe)G>,{e + c.) (34) 

Similarly, from the retarded and advanced components 
of the self-energy, the imaginary part is found to be 

''de 



and the lesser component 



Y^^L'/^[7'I<(£)G>,(£ + c.) (35) 



0,-yy' 2 F dS 

~ l67 ""7 2^'^ 



'^■I>{e)G<{e + uj). 



(36) 

Using the identities ^ EU it is now straightforward 
to carry out the integral over e, together with the sum- 
mation over a, a' and 7', to find that 

TT 

4 



Tf{~B/2) 



N- 



(2) 

n 



(37) 



and 



E<f (-i?/2) 



where the coefficients are given by 



(38) 



TV, 
TV, 



glj,V cotHV/2T) + [gl^ + gl^)T, (39) 



glMl + N{B)) 



(40) 



2(5iL 

+2glj,[{B + V)il + N{B + V)) 

+ {B-V){1 + N{B-V))], 
-2{glL + 9lR + 2glji)B, (41) 




'^^x^^ {uj) — J^qi j ^'iSiis) Q^i {e + U}) , (32) FIG. 3; Unlabeled Feynman diagrams for the p/ self-energies. 
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and N{lj) = l/(e"/'^ - 1) is the Bose-function. While 
the coefficients N^^i do not depend on A, i?*^^^ is in fact 
given by (n^ A — l)8g^i?. Since, however, the factor ni^\ 
vanishes in the Umit of A — > oo and leaves behind a term 
which remains finite, we can simply omit this term and 
employ the identity H4f |l for R. 



different orientations on the ce-loop. One involves two 
Cooper-bubbles and the other, two Peierls-bubbles, and 
the combined a6-component translates to 



2. Third order self- energy 



The third order self-energy diagram depicted in Fig.|3 
gives rise to two different terms, corresponding to two 



de 
2^ 



■d'c" 



cd' /~idc 



"d" a' v^d"c 



c'd /-idc 



(42) 



where a summation over ce-spin has produced the tensor 
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and O^^i^i' = —0^^,^,,. A contraction of the relevant Keldysh indices allows one to express as a sum of products 
of two vertex functions. Using again the relations derived in section FlI Al together with the analytical properties of 
I, the K component of the third order self-energy may be written as 



.K(3) 
^7A 



128i 



Re ^,I< {e)G>, ie + uj)+ {e)G<, {e + uj) 



iv 



-V), 



where we have introduced the shorthand J; 



3 

aa' a'' 



Jaa' Ja'a" Ja"a- Thc imaginary part takes the similar form 



.(3), 



_}_aP j3 

128 77'7" ckck'ck' 



and the lesser component reads 



Re [-^'lie)] (e)G> , (e + w) - {e)G<, (e + u) + {V 



-V). 



.<(3) 
"7A 



(^) — Y2g^^77'7"'^aa'a' 



-V). 



(43) 



Comparing to the second order result, the basic difference is the presence of the extra real part of the vertex-function, 
which provides the logarithmic enhancement underlying the Kondo effect. 

Again, the relevant integrations and summations are carried out using the specific form of the various components 
of I, and the result is expressed just as in Eqs. (|37(l and H38(l . but with new third order coefficients 



T +{V + B) coth 
D 



V 
2T 



In 



D 



+2{gtL + 9%R)T\n 



= ^9lRgd{ {l + N{B + V)) 



+ {l + N{B-V)) 
+ {l + N{B)) 
R^'^ = -8glRgdS^{B + V)\n 



[B + V) In 
{B - V) In 



\v 


+ B\ 




D 


\B 


+ v\ 




D 



T+{V - B) coth 



V 
2T 



In- 



D 



V~B 



\B + V\ 



Sin- 



D 

W\ 

D 



Vln 



D 

W\ 

D 



{B + V) In 
D 



D 



D 



B + V 



[B ~ V) In 



D 



D 



B-V 



B 



^iglL+g%R)B{l+N{B))\n 



D 



(44) 



(45) 
D 



^ In - 4(5iL + 9RR)Bil + NiB)) In --, 



B 



(46) 
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using the shorthand ln(D/|a;|) = ln(D/Va;2 + T^). As 
was the case for the second order terms, only i?^^^ de- 
pends on A, but in such a way that one can safely take 
the limit A — *■ oo which leads to H4()|) . 

Logarithmic correction of the form ln(Z?/|i?|), 
ln(D/|y|) or \n{D/\V ±B\), as in Eqs. appear 
throughout this paper. They are signatures of resonant 
scattering from one Fermi surface to another. More pre- 
cisely, a term like ln(_D/|F|) shows that resonant scatter- 
ing from the left, to the right Fermi surface is prohibited 
due to the difference V in their electro-chemical poten- 
tials. Similarly, a spin- flip process within the left lead is 
cut off by the required energy B, leading to h\{D/\B\). 
Most interesting are probably the logarithms of the form 
hi{D /\V — B\), characteristic for resonant scattering from 
the left, to the right lead, where the energy mismatch 
of the two Fermi surfaces is compensated by a spin-flip. 
These terms are especially important, insofar as they lead 
to pronounced cusps at ^ ~ i? in physical quantities (see 
below) . 



B. Solution of the Quantum Boltzmann Equation 

Having determined the pf self-energy in terms of the 
unknown nonequilihrium pf occupation numbers, the 
Keldysh-component of the Dyson equation provides a 
closed equation for n^\. Since we assume the system 
to be in a steady state, this equation is greatly simplified 
and may be expressed as^"^ 

r^;,H(7f,M=^^(c.)i]:f;,H, (47) 

or equivalently 

r^A(^)e<AM=AM5^<AM> (48) 

which merely states that the collision integral in the 
quantum Boltzmann equation has to vanish in a steady 
state situation. 

The pf spectral function appearing in this equation 
may be determined by solving the retarded and advanced 
components of the Dyson equation and takes the usual 
form 



(w) 



(w - Re[I]^(a))] -I- 7^/2)2 + {T^{uj)/2)- 



(49) 



where the shift Re[S] and the broadening F of the pf- 
energy levels are determined from perturbation theory, 
including leading logarithmic corrections. However, as- 
suming that A^{uj) is non-zero for all frequencies it can 
be divided out of Eq. which then takes the simple 
form 



(I]<,(c.) - T>^{uo))n^x{uj) = E<,(a;), 



or equivalently 



E>,(^)/S<,(c.) 



(50) 



(51) 
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FIG. 4: Nonequilihrium magnetization (dashed and long- 
dashed lines) and susceptibility (full and dotted lines) as func- 
tions of B/y, for T = 10" V and D = 10 V. Solid and 
dashed (dotted and long-dashed) lines correspond to 3.(2.) 
order perturbation theory. Upper panel: qll = Qrr ~ Qlh = 
0.05. Middle panel: qll = Qrr = 0.05 and qlr = 0.2. Lower 
panel: qll ~ Qrr ~ 0.2 and qlr ~ 0.05 (note the differ- 
ent scale for x)- Corrections from third order are more pro- 
nounced for gin / gLLQRR- 



In this equation the self-energies S^^^ are determined in 
perturbation theory as integrals involving the unknown 
distribution function n^\{uj). However, in bare perturba- 
tion theory, where there is no dressing of internal lines, 
n^\{uj) is always multiplied by the unrenormalized pf 
spectral function H12|l inside the integrals, and only the 
occupation numbers n-yx = ^7a(— 7^/2) will appear in 
S^'^. In our case, the quantum Boltzmann equation 
can thus be solved without any feedback from the re- 
tarded and advanced Dyson equations, and we can focus 
our attention on the on-shell occupation numbers n^x- 



Dividing out the spectral function and setting uj 
-jB/2, Eq. igHl) takes the form 



iT^xi-lB/2)n^x 



(-7S/2). 



(52) 



Since n^x, and thereby ^^x vanishes in the limit of 
A CX3, the whole equation may be divided by {Q)x and 
after this limit has been taken one can establish the same 
equation for the physical, projected, occupation numbers. 
Expressed in terms of the coefficients defined in Eqs. H37() 
and H38|l. we arrive at the equation 

iN^^+N^l+R)n^^N^^n^+Nuni, (53) 
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which can be viewed as a rate-equation with the transi- 

tion rates T^xT = f (^Ti + R) and Wj^ = fiV^x- This 
equation is readily solved together with the constraint 
equation n^- + = 1 and one finds that 



N- 



71 1 



n 



(54) 



or expressed in terms of the magnetization M = n| — 



Note that to obtain the observable magnetization to or- 
der g^\ii{D), it is not sufficient to consider only the on- 
shell occupations appearing in Eqs. (jSHSHl- In addition, 
one has to consider also contributions from the pf spec- 
tral function [see Eq. H14|l ] , which are discussed in detail 
in the appendix. 



M 



-R 



2N- 



n 



R 



(55) 



Inserting the second order expressions l|39H41fl , one ob- 
tains 



M{B, V) 



i9lL+9RR)Bcothi 



2TJ 



91r [{B + V) coth (^) + {B~V) coth (^)] 



(56) 



which, up to a factor of 2 in the definition of M, is exactly what was found in Eq. (4) of Ref. lSa. where a rate-equation 
like (|53(l was solved using second order transition rates. Including the third order corrections H44H46|I . we obtain 



with 



= '^9lR9d 



(y + B) In 



M ^ r/Af 



, ^ , -IV ~B) In , ^ , 
\V + B\ ^ ' \V-B\ 



B J2(9L + 9lJ ( 1 + 2ffaa In -|t 



\B\ 



coth 



-) 

2T 



'^9lR9d 



{V + B)hi- 



D 



coth 
' coth 



fV + B^ 




V 2T , 




/V-B^ 




V 2T , 


1 91r 



V + B\ 
(V + B) (^l + 2gd\n 

(V^B) f l + 2.9<jln 



{V - B) In ■ 
D 



V + B\ 
D 

V -B\ 

I 



D 


+ 


bE9L( 

Q. 


\^ 


V-B\\ 




2Vgd\ii 


D 

W\ 


+ 2Bgdln 


D 

\B\ 


2Vgd\ii 


D 

M 


-2Bgd\n 


D 

W\ 



D 

W\ 



(57) 



(58) 



(59) 



As expected, the Kondo effect reveals itself in logarith- 
mic enhancements, and interestingly enough the logarith- 
mic corrections to M come as [1 -I- (7ln(..)] rather than 
[1 + g'^ ln(..)], which is found in equilibrium. In the limit 
where T ^ V, the logarithmic corrections from T and JV 
cancel and we recover the usual thermal magnetization 
M = tanh.{B /2T). In the highly asymmetric case, where 
gLL 2> gRR,gLR, a similar cancellation takes place and 
one finds again that M = tanh(i3/2T). 

In the appendix, we demonstrate how the observ- 
able magnetization receives additional corrections of or- 
der g^ln(..), arising from self-energy corrections to the 
pf spectral function. Nevertheless, in the case where 
F 3> T, such additional corrections are subleading. 

In Fig. 21 we plot the magnetization and the cor- 
responding susceptibility as functions oi B/V for 
T ^ while Fig.jSjinvestigates the V dependence of M 
and dM/dV, which will both influence the conductance. 
Roughly speaking, the magnetization curve resembles the 
usual thermal magnetization with V replacing T, and the 



impurity-spin becomes polarized only when B exceeds V. 
However, structures close to ~ -B are much sharper 
and obtain quite different logarithmic corrections com- 
pared to the equilibrium case. The exchange-correlations 
change the slopes of the magnetization, which is clearly 
seen in the susceptibility (Fig.^ and in dM/dV (Fig.lSJ 
which become sharper spiked when including third order 
perturbations. In the limit of ^ max(i?, T), one may 
expand to find 

M^^lfds^) (l + 2g^^ln^) 

which simplifies to M = {2B/V){1 + 2g\n\V/B\) when 
gLL = 9rr = 9lr- Note that we still imply the cut-off by 
T meaning that \B\ and |V^| in the argument of the loga- 
rithm should be interpreted as \/B^ or + T^, 
respectively. 
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FIG. 5: Nonequilibrium magnetization as a function of V/B, 
for T = IQ-^B and D = lO^B. The two lower curves in 
each panel show the absolute value of the slope, i.e. |^^|. 
This quantity enters the conductance shown in Fig. |7| Solid 
(dotted) lines correspond to 3.(2.) order perturbation theory. 
Upper panel: qll = Qhr = Olr = 0.05. Middle panel: qll = 
Qrr = 0.05 and qlr = 0.2. Lower panel: qll = Qrr = 0.2 
and QLR = 0.05. 



It is interesting to note that the structure of loga- 
rithmic corrections is rather special in the case when 
the Hamiltonian ^ is derived from an underlying An- 
derson model by means of a Schrieffer-Wolf transforma- 
tion. In this case the exchange couplings are related as 
din ~ QLLgRR (cf. e.g. Ref. li^ and for such realizations 
of the model, H60|l simplifies to 



M«-(l + ^i^-t^ 

y I '2gLLgRB^ 



\V\' 



(61) 



disclosing a partial cancellation of logarithms, in the 
sense that the bandwidth, D, no longer enters the loga- 
rithmic correction [the same is also valid for the nonlinear 
magnetization (|57|) ]. This explains why the plots of the 
magnetization in the upper panels in Figs. 0] and [S] hardly 
differ from second to third order, while more pronounced 
effects are seen in the lower panels where g^j^ ^ gLLdRR- 

From Eq. H61() we see that the peak at _B = in the 
magnetic susceptibility, seen in Fig. 01 grows as ln(|T^|/T) 
for T ^V. Also, from 1)60(1 we learn that the high voltage 
tails in Fig. [S] fall off as 1/V. The derivative shown 
in Fig. |5l displays a peak near V — B already in second 
order perturbation theory which is then slightly enhanced 
or reduced by the third order correction, depending on 
the relative size of gd and gLR- As we shall see in section 
IIVI these features turn out to have a marked influence 
on the conductance. 

For arbitrary T and V , the susceptibility at _B = 
takes the following form: 



LR 



1 + 2gd In f + 2gd In ^) + ^„ g^^ + 2gL In f + 2.g2 ^ In ^ 



2V^coth (^) (1 + Again ^) + 2rE„5aa (saa + 2gL In f + 2<?i^,ln ^) 



(62) 



In the limit of T ^ V the logarithms all take the form 
ln(Z3/T); the corrections in numerator and denominator 
cancel, and we are left with the usual Curie law x = 
1/2T. In the extreme nonequilibrium situation, however, 
where V ^ T, the corrections no longer cancel and we 
are left with a complicated fraction times \/V. 

It is important to note that numerator and denomi- 
nator in H62|l have a rather different structure of loga- 
rithmic corrections, e.g. the g\^ term in the numerator 
receives corrections of the form 2(7d[ln(_D/T) -|-ln(I?/|l/|)] 
whereas the corresponding term in the denominator has 
the form Aqd h\{D l\V\). This observation was the basis 
of our claim in Ref. |4a that the perturbative renormal- 
ization group has to be formulated in terms of coupling 
functions which depend on the energy of the incoming 



electron. This will be explained in more detail in a forth- 
coming publication^. 

IV. NONLINEAR TUNNELING CURRENT 

The current operator measuring the charge flow from 
left to right lead is found from the equation of motion for 
the charge density in the left lead: 

dtUL ^ i[H,nL]. (63) 

The expectation value at time t is 

JL = iJLR^(^S{t) ■ [cl^,^,{t)f^'aCRka{t) - L R]j , 

k,k' 
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FIG. 6: Unlabeled Feynman diagrams for Dlr- Triangles 
denote the bare measurement vertices. 



and by defining a two-particle Keldysh contour-ordered 
correlation function 



DLRir,T') = 



(64) 



k,k' 



the current may be expressed simply by the imaginary 
part of the Keldysh component 



JL ^ JLRlm[Dfji{t,t)] 



(65) 



The lowest order contribution to this correlation func- 
tion is obtained from the first Feynman diagram in Fig.|Sl 
which translates to 



' > 8i 2tt 2tt 2tt 



+^ - 7. „ z::^ (66) 
y~ll,,.Qt^' {io)g'r-{Lo + f])A^^G^'=(e)Gf (e + 



with measuring vertices 7^ and 7^ defined in 1101 lll() 
and the spin-tensor from (jSSJ- Contracting the Keldysh 
indices and expressing things in terms of the vertex func- 
tion I, this simplifies to 



128z'^W^ 
X [^X^(£)^<X«(e)+^I^(£):«,I^(e)] , 

which in turn leads to 



.(2) ^ 



P!«/^Im[^I(£)]^'X^(£) (68) 

(y ^ -V). 



Evaluating this expression and performing the projec- 
tion, one finds the physical charge current to be 



iV + B)coth[^^] (69) 



-{V - B) coth 



2T 
V-B 
2T 



The magnetization entering this formula is the on-shell 
magnetization determined in section Fill Bl and since we 
include no dressing of internal lines this quantity receives 
no additional renormalization from the retarded and ad- 
vanced Dyson equations. 

The second order corrections to the correlation func- 
tion Dlh, are contained in the second diagram in Fig.El 
After contracting indices with the measurement vertices, 
this may be written in terms of Peierls and Cooper bub- 
bles as 



D 



K{2) 



LR '(^'*) — ^^^JLaJaR^^^^'f" j 

^'a^yK.ba'i^) J^ab" Y'KLb"a" (s) jK^'a (s) 



(70) 



AC'd Ll^dc ( \ \cd" a TP' d"c" f ^\ RT^c c / \ 

which upon full contraction leads to 



{^.T'^(e)^..X^(ef,X^{E) + ;^X^(e)^^X^(£)^Z^(e) 

R 

7'- 



(71) 



^■X^(£)^^,X^(£)^Z^(£)] 



(67) Inserting this into Eq. H65|l . we end up with 



.(3) _ 77'7" ^2 

•'^ 64i 



JIr{Jll + Jrr)J^ (72) 

Re [iflie) +^I{e)] Im [^^'lie)] i;!I^{e) 
-iV^~V), 



J 



which may finally be evaluated and projected to obtain the physical current 



JL 



TT 2 I../. . . D 



(V + B) l + 4.gdln. 



D 



V + B 



(F-B) l + 4.g<iln. 



D 



V-B\ 



(73) 



-M 



cothi^^^] l{V + B)(l + 2gd\n ^ 



2T 



coih\^^^ ) { {V-B) ( 1 + 2.9, In 



V + B\ 
D 

V -B\ 



2Fgdln^ + 2B5<ilnp 
-2Vgd\i^^-2Bgd\n^ 
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FIG. 7: Conductance as a function of V/B, for T = 10"^B 
and D = lO^B. Solid and dashed (dotted and long-dashed) 
lines correspond to 3.(2.) order perturbation theory. Dashed 
and long-dashed lines have M = 1 and correspond to Appel- 
baum's result, which neglects the V-dependence of M shown 
in Fig. 1^ Upper panel: qll = Qrr = Qlr = 0.05. Mid- 
dle panel: qll = Orr ~ 0.05 and qlr ~ 0.2. Lower panel: 
Qll = Qrr = 0.2 and qlr = 0.05. 



Here again, ln(£'/|x|) is a short-hand for 
hi{D / \fjP^+T^) as all logarithms are cut off by T. 
The current has acquired logarithmic corrections, which 
again enter as [1 -f (7ln(..)], and in zero magnetic field 
we recover the conductance obtained earher in Refs. |2^ 
and 23 (reinstalling e and K): 



G(F) = ^(|5L«)'3(l + 45,ln 



D 

\eV\ 



tanh(i3/2T), instead of the correct expression H57|) . and 
thereby left out a number of logarithmic corrections to 
the resulting conductance. These early results by Appel- 
baum are therefore only valid for a local moment cou- 
pled dominantly to one electrode, in which case qll 3> 
gRRTdLR and therefore indeed M — tanh(B/2T), as was 
pointed out in section Ull Bl Having solved the quantum 
Boltzmann equation in section Fill Bl including the lead- 
ing logarithmic corrections, we can now correct for this 
omission. In Fig. |7|we plot the conductance as a fimc- 
tion oiV/B, for T <C V , using both the nonequilibrium 
magnetization given by H57H59|I and the equilibrium value 
{M K, 1 since T <^ B), corresponding to Appelbaum's re- 
sult. In both cases a cusp develops at \V\ « \B\, beyond 
which spin-flip tunneling processes are energetically vi- 
able. Including the nonequilibrium magnetization, how- 
ever, leads to a pronounced enhancement of the cusps, 
coming from dM/dV (shown in Fig. |SJ, which is miss- 
ing in Appelbaum's result. In the lower panel, where 
9lr. < gd, Appelbaum's conductance is closer to the 
true curve, as for small gLR. one is closer to equilibrium 
idLR = is an equilibrium problem). 

Our result for the conductance appears to be similar to 
a plot by Sivan and Wingreen^ for the Anderson model, 
but since that work does not contain an explicit analyti- 
cal expression for the current and the magnetization it is 
difficult to compare the two results. Furthermore, loga- 
rithmic corrections to the quantum Boltzmann equation 
appear not to be included in their approach. 

In the high-temperature limit T ^ max(V,i?), there 
is no effect of V on the magnetization, and we recover 
Appelbaum's result for the current, or rather the re- 
sult he would have obtained by performing the final 
energy integrations in his expressions (39) and (60) of 
Ref. I23 The closed form expression presented in Ref. 
is the conductance rather than the current, and as an 
aside we note that differences from his formulas (where 
the T-dependence arises through hyperbolic tangents) to 
our Eq. H73|l (with hyperbolic cotangents) arise due to 
an (invalid) approximation by Appelbaum, when he ap- 
proximated the derivative of the Fermi function by a 5- 
function. Using the identity /(cj + e)(l — /(cj + O)) = 
N{e ~ fl)ifiuj + fl) ~ f{uj + e)) between Fermi, and 
Bose functions, we find coth(x) -I- xdx coth(a;) instead of 
tanh(x) giving rise to a different prefactor (2/3 instead 
of 1) for T larger than either V ^ B or \V — B\. 

In the limit of|F — i?| <C T<C the current may 
be expanded to give 



JL 



gi^B l + 45dln 



D 
2LBI 



(75) 



(74) 



which remains finite as T — + 0. When expanded in the 
same limit, the conductance, however, takes the form 



In a finite magnetic field, however, the proliferation 
of logarithms is more intricate. While the logarithmic 
structure in H73|) was derived already by AppelbaumSSiS^, 
he tacitly employed the equilibrium magnetization M = 



G{B,T)^-gl^ 

4'' 



25d 



In - — - 

LB 



ln^-1 
T 



(76) 
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which reveals a logarithmic divergence as T ^ 0. 

This divergence seems to suggest that logarithmic cor- 
rections may grow large, signaling a transition to the 
strong coupling regime. Such behavior is, however, pro- 
hibited by spin relaxation processes emanating from the 
finite current flowing through the dot. As discussed in 
Refs. Is^ and li^ the relevant cutoff energy is the voltage 
dependent transverse spin relaxation rate At tem- 
peratures T < T2 the factor T in ln(_D/r) is replaced 
by r2 ~ g^V, leading to a finite correction even in the 
limit T ^ 0. This cut-off arises as a joint effect of vertex 
and self-energy corrections, which will be worked out in 
detail in a forthcoming publication'^-. 

V. DISCUSSION 

In the present work we have calculated the local mag- 
netization and the charge current through a quantum dot 
at large bias voltage and in the presence of a magnetic 
field. We have considered a Kondo model where the dot 
is represented by a quantum spin S — 1/2, coupled to 
leads by exchange interaction and tunneling. 

Remarkably, the structure of perturbation theory in 
steady-state nonequilibrium is rather different to that in 
equilibrium. The main physical reason is that in our 
problem the occupation on the dot is completely unde- 
termined in the limit of vanishing couplings to the leads 
and therefore has to be calculated from the solution of the 
quantum Boltzmann equation (i.e. from a self-consistent 
Dyson equation). More generally, in equilibrium all dis- 
tribution functions are given exactly by the "bare" Boltz- 
mann, Fermi or Bose functions without any corrections 
from interactions while out of equilibrium these functions 
have to be calculated and will depend on all the cou- 
plings. Probably the most drastic consequence is that 
for finite voltages the magnetization H56() is modified even 
for vanishing couplings, i.e., in zeroth order of perturba- 
tion theory as has been emphasized in Ref. |5^ As a 
consequence the structure of logarithmic corrections to 
the magnetization is also rather different out of equilib- 
rium. As the matrix elements of order in the quantum 
Boltzmann equation get corrections of order ln(D), the 
perturbative magnetization (|57() is generically of the form 

g2^g3 ~ const. + 0[gln(..)], to be compared to the 
equilibrium case where logarithmic corrections arise only 
to order g^ ln(..). 

In this paper, we have calculated the magnetization 
and the current through a Kondo dot to third order in 
the coupling, including the leading logarithmic terms. In 
contrast to earlier treatments of this problemSSiS^, the ef- 
fect of the nonequilibrium magnetization on the current 
is incorporated and indeed shown to be important for a 
typical quantum-dot experiment in the Kondo regime. In 
finite magnetic field the differential conductance exhibits 
threshold behavior &i V = ±i3, reflecting the fact that 
spin-flip tunneling is possible only for \V\ > \B\. Tak- 
ing the voltage dependence of the magnetization properly 



into account, the conductance shows a cusp at V = ±_B, 
already to lowest order in perturbation theory. Going 
one order higher, logarithmic corrections are found to 
enhance these cusps even further, increasing the conduc- 
tance substantially over the threshold plateau, even for 
magnetic flelds much larger than temperature. This be- 
havior clearly calls for a resummation of the perturbation 
series to inflnite order. 

The intricate structure of logarithmic corrections, re- 
vealed by nonequilibrium perturbation theory, enforces 
a modified formulation of the perturbative renormaliza- 
tion group. Thinking in terms of poor man's scaling, a 
resummation of the logarithmic corrections cannot sim- 
ply be achieved by collecting the logarithmic corrections 
into a renormalized coupling constant. This is clearly il- 
lustrated by Eq. (|62|l for the susceptibility, which shows 
that different combinations of logarithms appear in the 
numerator and in the denominator. However, as we have 
recently demonstrated^, exactly this structure of leading 
logarithms is generated in a poor man's scaling approach 
dealing with energy dependent coupling functions. A de- 
tailed account of this approach will be given in a subse- 
quent publication^. 

An additional important difference between nonequi- 
librium and equilibrium physics is the fact that in a 
nonequilibrium situation the current through the system 
generates substantial noise. As a consequence, quantum 
coherence is limited to energies above a certain relaxation 
rate T2. The scale r2 will therefore cut off all logarith- 
mic divergences which remain in the limit T —^ in the 
perturbative expressions for magnetization and current 
found in this paper. This important piece of physics is 
not included in the low-order bare perturbation theory 
presented here and requires a resummation of subleading 
self-energy and vertex corrections. This will be demon- 
strated explicitly in a forthcoming publication^. 

Historically, the Kondo effect has played an impor- 
tant role in the development of techniques like the renor- 
malization group^'^ to treat strong-coupling problems for 
systems in thermal equilibrium. It is our hope that the 
present perturbative calculation of leading logarithmic 
corrections can serve both as a starting point and a check 
for future developments of similar methods applicable to 
systems out of equilibrium. 
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APPENDIX: MAGNETIZATION AND 
SELF-ENERGY CORRECTIONS 



numbers. Inserting (|A.4|I into Eq. (|51|l shows that indeed 



In this appendix we calculate the observable magneti- 
zation to order (7^1n(..), by including self-energy correc- 
tions to the pseudo fermion spectral function. Using the 
prescription for evaluation of a canonical ensemble 
average, the magnetization is determined as 



M= hm 

A^oo n|A + nix 



where 



T-7A 



(A.l) 



(A.2) 



and in principle, the full frequency dependence of both 
the spectral, and the distribution function is needed to 
carry out this integral. 

To illustrate how such renormalization works out 
within the pseudo fermion approach, we shall commence 
with the simpler case of thermal equilibrium and con- 
sider merely linear response for B <^ T. Traditionally, 
most work regarding the Kondo effect on the magnetic 
susceptibility in this regime has been conducted using a 
Kubo-formulaSiSSiSi, and to the best of our knowledge 
only one work has taken the approach outlined above to 
calculate the magnetization directly from the renormal- 
ization of the pf spectral function^. Since, however, this 
latter work contains an error, and as we think that the 
calculation is instructive, we shall consider this simpler 
case in some detail before briefly discussing the nonequi- 
librium case. 



1. Equilibrium Magnetization 

At zero bias-voltage the local spin is in equilibrium 
with the leads, and the pf distribution function reduces 
to a simple Fermi function. In the Keldysh formal- 
ism, this may be viewed as a simple consequence of the 
KMS (Kubo-Martin-Schwinger) boundary condition^—, 
which states that in thermal equilibrium Q^-^{uj) = 
—Q'^-^{uj) exp[(a; + A)/r]. The quantum Boltzmann equa- 
tion igHI), may be rewritten as 



^^{Lo)g<^{u) = ^<^{Lo)g>^{uj), (A.3) 

from which the KMS condition is seen to imply that 



(A.4) 



This can also be verified explicitly from our perturbative 
expressions for S^'*^, by employing the (KMS) condition 
1 — n-yA = n-yAexp[(A — ^B/2)/T] (exphcitly satisfied by 
a Fermi function) for the on-shell equilibrium occupation 



(A.5) 



and when derived in this way, it becomes clear that the 
KMS condition ensures a highly non-trivial cancellation 
of interaction corrections in the case of thermal equilib- 
rium. Applying a finite voltage, this condition is vio- 
lated and the nonequilibrium distribution function will 
be affected by interactions in the manner which we have 
described in section ITlI Bl 

As mentioned earlier, the first order pf Hartree self- 
energy vanishes unless one includes a Zeeman term for 
the conduction electrons. However, adding such a term. 



HI 



Zeeman 



B 



E 

T;k' cj' 



,C / 



(A.6) 



to the Hamiltonian, one finds that the self-energy is en- 
tirely real and given as 



= lBgd/2, 



(A.7) 



where the ce Green-function now depends on the spin. 
Note that the Keldysh component of the Hartree self- 
energy is identically zero, due to the fact that Re [G'Q.y(e)] 
is an uneven function. In a case where particle-hole sym- 
metry is broken on an energy-scale S, one obtains a finite 
contribution of order J5/D, which can be neglected for 
large D. The only effect of including a Zeeman term for 
the conduction electrons is therefore a constant shift of 
the pf energy levels. 

The second order pf self-energy has a real part, which 
has only a negligible contribution from the term 
Without this term, one finds that 



Re[I](2)(c.)] = - 



0^ 

■Jaa 

32 

X In 



D 

de tanh 

-D 



e + UJ — ^a' 



2T 



which implies that 



and 



Repf (-7S/2)] 



a^Re[I](2)(c^)] 



9 D 



3 2, D 



(A., 



(A.9) 



(A.IO) 



in the case of F = and T B. 

The broadening of the pf energy levels is given by the 
imaginary part of the self-energy asr = i(S]> — 1]<), and 
when neglecting E^, which projects to zero, we end up 
with r — iS^ which is given by 

r^(^) = g^e^yiiu + 7'S/2)(1 + N{lu + 7'S/2)), (A.ll) 
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to second order in g. This function is highly asymmet- 
ric and expanding to first order in B/T one finds the 
asymptotic behavior 



to which one may finally add the corresponding induced 
spin-polarization of the conduction electrons to obtain 
the total magnetization 




7_B 
6T 



(A.12) 

complemented by the fact that r^(w) = for > 2D, 
since the excitation of particle-hole pairs giving rise to 
the broadening is limited by the bandwidth. Since F is 
essentially constant near the peak of the spectral func- 
tion, one may approximate Aj in (|49|l by a Lorentzian for 
LU € [-T,T], and by T^{uj)/uj'^ for \uj\ > T. Introducing 
the wave function renormalization factor 



= 1 Q f n — , 

2 T' 



(A.13) 



the coherent part of this approximate spectral function 
integrates to Z on the interval [—T,T]. The exponential 
integral of the negative frequency tail contributes with a 
number of the order of g^, which should be neglected, 
while the integral from T to D yields exactly 1 — Z, 
which ensures that this approximate spectral function 
integrates to 1. 

To find the magnetization, the integral (|A.2p may now 
be evaluated using the Boltzmann distribution and this 
approximate spectral function. The spectral function is 
centered at a frequency coj satisfying the equation ujj = 
—jB/2 + Re[T,j{ujj)], and consequently the integral over 
[— T, T] contributes a factor of Zn~^\{uj~^'), that is 



27r - uj-^Y -I- (r^/2) 
1-.9 



^Zn^\{uj) 



Z + 



7 2, D 
- a in — 
2 .y J. 



(A.14) 



to first order in B/T. Since the distribution function falls 
off exponentially for uj ^ T, the integral over [T, D] is 
negligible. On the negative frequency tail, the spectral 
function decays exponentially, but this is compensated in 
the integral by the exponential increase of the Boltzmann 
distribution, and one finds that 



duj r. 



2, ^ 



(A.15) 



Adding up (jA.14p and (|A.15|) and inserting in l|A.l|) . one 

ends up with the magnetization 



M 



2T \ ^ ^ T 



(A.16) 



(5^ 



\ Pauli 



^("1-29-492 In 



(A.17) 



Here the z-component of the total spin of lead a has been 
introduced as 



k,cr;k' ,cr' ;a 



(A.18) 



and the pf and ce g-factors are assumed to be equal. 
This result matches the high-temperature expansion of 
the exact Bethe ansatz solution as it shouldiiSS. 

In Ref. I52L the normalization of the approximate spec- 
tral function is demonstrated just like here. In calculat- 
ing the magnetization, however, the contribution IjA.lSp 
was not included and the Z-factor was argued to be can- 
celed by the same factor appearing in the denominator 
of ljA.l|) . Altogether, this error leads to a prefactor of 2 
instead of 4 in front of the 5^ ln(£)/r)-term in (|A.16p . 
which destroys the correspondence with the exact result. 



2. Nonequilibrium Magnetization 

As demonstrated in the equilibrium case, the magneti- 
zation is renormalized in a delicate balance between shifts 
and broadening of the pf energy levels, i.e. between the 
influence of Rc[S]] on the coherent part, and of F on the 
incoherent tails of the spectral function. In the case of 
finite bias-voltage, however, the renormalization of the 
distribution function becomes important. 

At finite voltage, there are logarithmic corrections to 
and S^, which no longer cancel in Eq. H51|) . and 
the magnetization now exhibits the stronger renormal- 
ization by factors of g\n{D) rather than g and g^ \n{D). 
To properly describe the crossover to equilibrium, as T 
becomes greater than V, one should include these sub- 
leading corrections deriving from the renormalization of 
the pf spectral function. This can be done in much the 
same way as above, as long as care is taken to separate the 
coherent part of the spectral function from the incoher- 
ent tails at either T or V, so as to ensure normalization 
of the total spectral weight. When T is increased beyond 
V, the .gln(_D) corrections get less and less important. 

It should be emphasized that, as mentioned already in 
section Hvl the effects of shifts and broadening, discussed 
in this appendix, affect only the observable magnetiza- 
tion and have no influence on our result for the leading 
logarithmic corrections to the current. 
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